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We have investigated Haldane's conjecture for the S=2 isotropic antiferromag- 
netic quantum spin chain with nearest-neighbor exchange J. Using a density 
matrix renormalization group algorithm for chains up to L=350 spins, we find in 
the thermodynamic limit a finite gap of A = 0.085(5)J and a finite spin-spin cor- 
relation length £ = 49(1) lattice spacings. We establish the ground state energy 
per bond to be E = — 4.761248(1)J. We show that the ground state has a hidden 
topological order that is revealed in a nonlocal string correlation function. This 
means that the physics of the S=2 chain can be captured by a valence-bond solid 
description. We also observe effective free spin-1 states at the ends of an open 
S=2 chain. 
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Since Haldane's conjecture |1| that the physical properties of (isotropic) anti- 
ferromagnetic quantum spin chains depend crucially on whether the spin is integer 
or half-integer, his prediction of a ground state with a finite spin-spin correla- 
tion length and a finite gap to spin excitations for integer spins has been studied 
in numerous papers. Numerical methods have served to establish quantitative 
results where analytical work had to rely on often uncontrolled approximations. 

Both the existence of a finite gap between a ground state singlet and excited 
states and of a finite correlation length is well established for the S=l isotropic an- 
tiferromagnetic quantum spin chain|2[]. Recent estimates of the Haldane gap are 
discussing the fifth digit: A = 0.41050(2)0 from the Density Matrix Renormal- 
ization Group (DMRG) algorithm[§J, and A = 0.41049(2) f| from extrapolation 
of exact diagonalization results for chains up to length L=22. Similarily, the 
correlation length is £ = 6.03(1) || ||. The spin wave velocity is then found to 
be c = A£ = 2.475(5), to be compared to the semiclassical value c = 2S = 2. 

In this work we have established numerical values for the respective quantities 
for the S=2 isotropic antiferromagnetic spin chain with nearest-neighbor exchange 
using the DMRG algorithm. We show that Haldane's conjecture is obeyed in this 
case: the Haldane gap also exists for S=2, further supporting the case that the 
gap exists for all integer spins. In addition our results show that a valence-bond- 
solid picture^ is valid for the S=2 spin chain. 

The investigation of the S = 2 case is hindered by the fact that one expects 
As =2 <C As = i and £s =2 3> £5=1. This belief is due to the asymptotic form A s ~ 
S 2 exp(— 7^) [[!]]. One may approximatively expect A s=2 ~ 0.07 and £s=2 ~ 70. 
The much larger correlation length makes finite size extrapolations feasible only 
for much longer chains than in the S=l case; furthermore, the inherent statistical 
or systematical imprecisions of all numerical methods become more worrisome 
due to the small size of the gap. At the same time, the number of states per site 
rises from 3 to 5, greatly reducing the length of numerically tractable chains. At 
the present time, there are only a few preliminary estimates ||, |§ in the literature. 

The DMRG algorithm allows to treat very long chains while retaining good 
precision for the energies and expectation values by an iterative truncation of the 
Hilbert space down to M basis states. For a detailed discussion of the method 
we refer to |3], ||. In the following we have investigated chains up to a length 
of L=350, which we will find to be about 7£. This allows for sensible finite size 
extrapolations. At the same time, we have retained up to M=210 states, which 
allows for reasonable accuracy. As the DMRG is much more precise for open 
than for closed chains, we have considered only open chains. 

To calculate the gap, we consider an open S=2 chain with a spin-1 at each 
end (Si — Sl — 1; Si — 2 otherwise): 

L-2 

H = Jend Sl • S2 + J Si ■ Sj+i + J en d Sl~1 ' Sl- (1) 

i=2 

We set Jend = J and the results will be given in units of J from now on. 
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Throughout this study we are guided by the so-called valence-bond-solid 
(VBS) wavefunctionsj?]]. They are exact ground states of very special spin Hamil- 
tonians (VBS hereafter). For S=l it is well established that the simplest Heisen- 
berg Hamiltonian fl(| shares a similar physics with the VBS Hamiltonian. These 
VBS ground states have a finite correlation length and there is a gap to spin ex- 
citations. From these VBS Hamiltonians it is known that the ground state of the 
open VBS-chain with integer spin 5 is exactly (5+l) 2 -fold degenerate 0, because 
the chain ends act like free 5/2 spins. For a generic "Haldane-like" Hamiltonian 
we expect this manifold of states to be split due to interactions between end spins. 
We couple a spin-5/2 to each end with an antiferromagnetic coupling J en d > 
(we chose J en d = 1) to select a state in this manifold. 

We calculated the gap for L=60,90, 120, 150,210,270, as the difference between 
the total energies of the lowest energy states in the S^ otal = and S^ otal = 1 
subspaces. We have kept M=210 states in the DMRG. White has given a 
procedure^, which allows to ameliorate iteratively the results found for the to- 
tal energies of a chain of fixed length L. For short chains, the DMRG algorithm 
truncates less well than for long chains; after the chain has grown to full length 
L, the results can be used to ameliorate the truncations at the start of the algo- 
rithm. We do not apply this iteration to the ground state, but iterate once for 
the excited state. From test calculations we found that iterating (further) does 
not substantially increase the precision obtained or influence our extrapolation. 
We extrapolate the results first for fixed L to M=oo, and then these extrapolated 
values are used to investigate the thermodynamic limit L=oo. 

We have extrapolated the energies to M=oo and deduced the gap values by 
subtraction. We find that both for the (uniterated) ground state energy and (it- 
erated) excitation energy the deviation from the asymptotic value is proportional 
to (1 — J2 Aj), where A, are the eigenvalues of the states of the density matrix that 
are kept at the last iteration. The proportionality coefficient is in both cases of 
order 10 (much larger for the uniterated excitation energy), the truncation error 
1 — Aj for M=210 between 3.7 x 10~ 8 and 9.4 x 10~ 8 for the ground state energy 
and 6.7 x 10~ 8 and 5.3 x 10 -7 for the iterated excitation energy for L between 
60 and 210. The precision of the gap is in reality better than estimated, because 
both energies are overestimated due to the variational character of the DMRG, 
and there will be partial cancellation of the errors. 

The ground state energy results allow us to calculate the ground state energy 
per site E . We find that for L«120, for each M< 210 considered the ground 
state energy per site has converged to a precision of 10~ 6 . These values we then 
extrapolate in M, again finding a linear dependence of the error on the truncation 
error. We conclude that E = -4.761248(1). 

We can extrapolate our results for 90 < L < 270 to L=oo by fitting an L _1 -law 
(figure 1) and obtain a gap estimate of 0.081(1) for L=oo. We note that deviations 
from perfect linear behaviour are extremely small in the range of chain length we 
have considered. However it is important to note that the asymptotic behaviour 
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of the gap in an open chain is expected to be 1/L 2 : the massive quasiparticles 
have dispersion 

E k = ^A 2 ^ + c 2 (k - 7r) 2 pa Aoo + c 2 (k - 7r) 2 /2A 0O near the bottom of the 
band. On an open chain the quasiparticle cannot stay in a stationary state 
with k = 7T but instead k — n pa ty/L as a particle in a box since there is no 
translational invariance. This means that the lowest excited state has a gap 
+ c 2 7r 2 /2A 0O L 2 + 0(1/ L 3 ). This behaviour has not yet been reached for 
L=270. Thus there must be a crossover point to parabolic behaviour of A(L). 
As a consequence we can obtain an upper bound on the gap by assuming that the 
parabolic behaviour sets in immediately beyond L=270. We match a parabolic 
curve A(L) = A^ + aL~ 2 to our extrapolated linear gap curve at L=270 such 
that up to the first derivative the two regimes meet continuously. This leads 
to Aoo — 0.090 (and a ^620). We have also estimated the correlation length 
(see below) to be ^50 and thus c = £Aqo pa 4.2. This gives an estimate of the 
coefficient of the 1/L 2 parabolic term. Using this estimate (do ~1024) we now 
say that there is a crossover length L at which asymptotic behaviour A(L) = 
Aqo + clo/L 2 sets in. Matching the parabolic curve to the straight line requires 
Aoo = 0.085 and Lq pa 450. This is a perfectly consistent set of results. We thus 
choose as a central value for the gap Aqo = 0.085 the lower bound being fixed 
by the linear fit of our data and the upper bound being 0.09. We quote our final 
result as A = 0.085(5). A gapless state is excluded beyond any doubt. 

To calculate the correlation length, we partially lift the ground state degen- 
eracy of the open chain by adding a spin-1 on one end of the S=2 chain (Sl = 1; 
Si = 2 otherwise): 

L-2 

H = J Si ■ Sj+i + JendSL-l • Sl- (2) 

1=1 

The ground state will then be a spin-1 triplet; at the spin-2 end of the chain 
(position 1) one expects an effective free spin-1. If we consider the ground state 
with S^ otal — 1, it is expected that (S?) decays purely exponentially from a value 
around ±1: (Sf) oc (— l) 1 ^ 1 exp(— (i — l)/£)- From this equation the correlation 
length £ can be obtained much more precisely than from (Sj • Sj), as it is de- 
rived from a one-point correlation, where the systematic errors of the DMRG 
tend to cancel partially, whereas they build up for the two-point correlations, as 
correlated states in each half-chain are systematically neglected at the same time, 
underestimating £. 

Because of the slow increase of precision with M for S=2, we have to perform 
an extrapolation of £ with respect to M. At the same time we consider chains of 
length L=270 and repeat some calculations for L=350, to study possible finite 
size effects. We find that those are minor, as are effects of the choice of J enc i, 
and concentrate on the extrapolation in M. In figure 2, we give In \{S*)\ vs. %. 
The generic behaviour is the same for all M: After the expected decay of S z , 



3 



there is a small increase towards the other end of the chain (greatly exaggerated 
by the logarithmic scale). With increasing M, the minimum shifts to the right, 
the spin expectation value at the right end of the chain is greatly reduced, the 
decay becomes more exponential and approaches a limiting curve which we take 
to be the M=oo result. To derive £, we give in figure 3 the local decay length 
2/(ln|S'f| — In 1 S* +2 1) (averaged over two neighbouring sites to reduce odd-even 
site oscillations). One sees that at the spin-2 end of the chain, the correlation 
length is first rather short, but saturates to its bulk value (the shoulder). For 
M=180, the convergence of £ for M— > oo can be well established, and we find a 
correlation length £ = 49(1) in a system which has thus length L rs 5.5£, which 
is consistent with minor finite size effects. At the same time, we find that (Sf) 
converges towards 1.13(1), supporting the picture of an effective free spin-1 at 
the end of the chain. 

Furthermore, we have measured the following nonlocal string correlation: 



_ m— 1 

G(n,m) = (^exp(i- £ Sf)^). (3) 

l 



For the S=2 VBS-chain G(n,m) — > —1 for |n — m\ — > oo, revealing a hidden 
topological order. Investigating isotropic Heisenberg S=2 chains up to L=270 
and keeping up to M=150 states, we find that G(n, m) saturates to —0.726(2) for 
\n — m\ — > oo indicative of hidden topological long-range order in the Heisenberg 
S=2 chain|Tl|]. These results imply that the physics of the Heisenberg S=2 chain 
can be captured by a VBS description. 

For the isotropic AFM quantum spin-2 chain we find from A = 0.085(5) and 
£ = 49(1) a spin wave velocity c = 4.2(3). This result indicates that, as expected, 
S=2 is more classical than S=l. We have also shown that the S=2 spin chain has 
hidden long-range order as in the VBS model. There is excellent evidence that 



the S=l Haldane gap is realized in the compound NENP p2| . It is an open and 
interesting question to find an experimental system with S=2 spins to check our 
findings. 

It is a pleasure to thank N. Elstner, O. Golinelli and R. Lacaze for valuable 
discussions. The calculations were carried out on the CRAY C94 of the CEA. 



1 Figure Captions 

Figure 1: Gap vs. L' 1 for L=90, L=120, L=150, L=210 and L=270 for different 
M, the number of kept states. Full squares: extrapolated values for M— >• oo. Solid 
line: extrapolation of the gap values to oo, yielding a lower bound. The 
dashed line is an upper bound on gap behaviour which is obtained by assuming 
that the asymptotic form A(L) = A^ + a/L 2 starts at L=270. Since we have no 
evidence for such a behaviour, the dashed curve should be a strict upper bound 
on the gap. 
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Figure 2: In \ (S?) \ vs. site % for an open S=2 chain with a spin-1 attached at 
the right end (J en d — !)■ Note that for increased chain length the straight-line 
fit for the slope is possible for larger i, while the slope itself remains virtually 
unchanged. The expectation values for the right end of the chain have obviously 
not converged yet, whereas the asymptotic oo slope value can be determined. 

Figure 3: Two-site average of the local correlation length 2/(ln | (Sf) \ — 
\n.\(S? +2 )\) vs. site i for i < 125 for the same chain as in Figure 2. The de- 
viations from the shoulder in the center of the diagram are finite size and finite 
precision effects (for the right end). 
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